Rupture of a liposomal vesicle 
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We discuss pore dynamics in osmotically stressed vesicles. A set of equations which govern the 
liposomal size, internal solute concentration, and pore diameter is solved numerically. We find 
that dependent on the internal solute concentration and vesicle size, liposomes can stay pore-free, 
nucleate a short lived pore, or nucleate a long-lived pore. The phase diagram of pore stability is 
constructed, and the different scaling regimes are deduced analytically. 



I. INTRODUCTION 

Liposomal vesicle consists of a lipid bilayer separating 
the interior volume, containing an aqueous solution, from 
an exterior suspension. The vesicle membrane allows for 
a free exchange of water between the interior and the 
exterior of the liposome, with the flux determined by 
the membrane composition. On the other hand, lipidic 
membrane strongly inhibits passage of large molecules, 
in particular, if they contain ionized groups. 

Liposomes are of great theoretical interest as the sim- 
plest model of a biological cell. They are also of great 
practical importance as vehicles for drug delivery. In the 
latter case liposomes are designed to contain a specific 
drug or a gene needed to fight the decease. The liposomal 
affinity for infected tissue can be increased by varying the 
membrane composition or including ligands which bind 
to specific receptors. 

If a vesicle containing high internal solute concentra- 
tion is placed inside a dilute solution, the osmotic influx 
of solvent into the interior of a vesicle can lead to its 
rupture. Whether, the rupture occurs depends on the 
membrane elasticity and on the internal solute concen- 
tration of the liposome. Rupture of the liposomal mem- 
brane results in formation of pores . This releases the 
membrane stress, but comes at a price of exposing the 
hydrophobic membrane interior (lipidic tails) to water. 
Once a pore is formed, the internal content of the vesi- 
cle begins to leak out, resulting in a decrease of mem- 
brane tension and eventual pore closure. We find that 
depending on the vesicle size and internal concentration 
of solute, pores can be either short or long lived. For 
long lived pores a scaling relation between the life-time 
of a pore and the size of the vesicle is found. The full 
phase diagram of pore stability in the concentration - 
vesicle size plane is constructed, and the different scaling 
regimes are deduced analytically. 

The paper is organized as follows: In section [H] we 
review, the previously derived equations governing the 
nucleation and growth of a pore in an osmotically stressed 
vesicle 0. In section Hill numerical solution of dynamical 
equations is presented. In section llVl the phase diagram 
for different dynamical regimes is derived. In section Ivl 
the rate of solute leak-out is determined and analytical 
estimates of the pore life-time are provided. Finally, in 
section lv*Tl the conclusions are presented. 



II. 



THE MODEL 



As was already stressed in the introduction, the li- 
posomal membrane allows for a free exchange of water 
between the exterior and the interior of a vesicle. The 
rate of this exchange is determined by the permeability 
of membrane, P. On the other hand, lipidic membrane 
strongly inhibits exchange of solute molecules between 
the inside and the outside of a liposome. When a vesicle 
of high internal solute concentration is placed inside a 
solute-depleted medium, an osmotic pressure difference 
causes an influx of water into the vesicle. A vesicle then 
swells until the internal Laplace pressure is able to com- 
pensate the the osmotic pressure. The influx of water 
result in a build up of membrane stress which energet- 
ically favors membrane rupture and formation of pores. 
Pores are nucleated in the membrane through thermal 
fluctuations. Here we consider the opening of a single 
pore. The underlying assumption is that, once a pore is 
formed, stress is quickly released and the creation of a 
second pore becomes highly unlikely. This situation is 
quite different from what is encountered in electropora- 
tion. In that case opening of a pore does not fully release 
the membrane stress, which is induced by the transmem- 
brane potentiaLand one finds a coexistence of pores with 
different sizes Jl2? ]. 

The single pore assumption allow us to write simple 
equations governing the internal vesicle dynamics. Des- 
ignating the difference between the internal and the ex- 
ternal molar concentrations of solute as c - and consid- 
ering, for mathematical simplicity, a spherical vesicle of 
radius R, and a circular pore of radius r — the mass 
conservation leads to 



4irp R' 



dR 
~dt 



j w - nr pv , 



(1) 



where p is the density of water, j w the osmotic current 
and v is the leak-out velocity. 

The osmotic current j w is determined by the perme- 
ability of the liposomal membrane and the difference be- 
tween the target osmotic pressure 



A Po = k B TN A c 
and the Laplace pressure, 
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inside and outside the vesicle. In the above expressions 
kg is the Boltzmann constant, T is temperature, Na is 
the Avogadro number, and a is the membrane surface 
tension. A simple phcnomcnological expression for the 
osmotic current of water into the vesicle is, 



j w = P(AnR 2 - 7rr 2 ) 



A Pl 



10 3 k B TN A 



(4) 



where the conversion factor 10 3 accounts for the use of 
molar concentration of solute c. 

If c is not too large, the membrane integrity will not be 
compromised, and a stationary state with j w = will be 
achieved. Under these conditions the osmotic pressure 
is completely compensated by the Laplace pressure, re- 
sulting in a zero net flux of solvent. For sufficiently large 
internal concentration of solute, a stationary state will 
not be achieved before membrane ruptures. The leak-out 
velocity HSH® 

of the internal content of a liposome is 
determined by the balance between the shear stress, pro- 
portional to rjv/r, and the Laplace pressure inside the 
vesicle ApL- For low Reynolds numbers Q 



Ap L r 
3irr] 



(5) 



where r\ is the solvent viscosity. 

The growth of a pore is controlled by the rate at which 
the membrane elastic energy is dissipated. Since the vis- 
cosity of membrane is five orders of magnitude larger 
than that of water, most of the energy dissipation is con- 
fined to the membrane interior Q, 



dr _ dE 
Vm dt ~ ~~dr~' 



(0) 



where I is the membrane width and r\ m is the membrane 
viscosity. A lipid bilayer has low permeability to solute 
particles, in particular if they are charged, so that the 
internal solute concentration is modified only through the 
osmotic influx of solvent or the efflux of solute through 
an open pore, after the membrane has ruptured. The 
continuity equation expressing this is 



3 dt 



AnR 2 c 



dR 
~dt 



irr 2 cv 



(7) 



where we have assumed that solute is uniformly dis- 
tributed inside the vesicle. In the absence of a pore, 
efflux is zero, and the second term on the right hand side 
of Eqs. (JIJ and Q disappears. 



A. The membrane energy 

The membrane energy consists of two terms. The elas- 
tic term E s , measuring the cost of increasing the mem- 
brane area beyond its equilibrium unstretched size Aq, 
and the pore contribution E p resulting from the partial 
exposure of the hydrophobic lipidic tales to the aqueous 
environment. 



Parameter Value 



Source 



7 

K 

f/m 
Vw 
P 



1(T 12 J/m Ref. [6] 

0.2 J/m 2 Ref. [9] 

100 Pas Ref. [6] 

0.001 Pas 

1.8 x 10~ 4 kg/(m 2 sM) Ref. [? ; 
3.5 nm Ref. \Q 



TABLE I: Characteristic values for the physical parameteres 
used in the calculations. 



For large osmotic pressures, which are of interest to us, 
the membrane thermal undulations can be ignored and 
the membrane elastic energy takes a Hooke-like form 



E s (R,r) = ^-n(A-A Q ) 2 



(8) 



where Aq = AttRq is the equilibrium surface area of 



an unstretched vesicle, A 



AttR 2 



TIT 



2 is the total 



membrane area, and k is the membrane elastic modu- 
lus H, 13, [13 . The membrane surface tension is 



dE s 
a = -dA 

and the pore energy is 

E P (r) 



A-A 

An 



2 7r 7 r , 



(9) 



(10) 



where 7 is the pore line tension. 

The typical values for the physical constants involved 
in the model are given in the Table 1. 



B. The rupture condition 

The total energy of a membrane containing a pore is 
E(R, r) — E S (R, r) + E p (r). A cost of opening a pore of 
radius r is then, 



AE(R, r) = E(R, r) - E(R, 0) 



(11) 



In Fig. 2] we plot AE(R, r) as a function of r for various 
ratios of R/Rq. For R ~ Ro, the membrane is relaxed 
and r — is the only minimum of AE. For R bigger than 
the critical radius R Cl the energy cost function develops 
a barrier located at 



n 



Rq 



v/3 

and a new minimum at 
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where 



ip(R, Rq) = cos 
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The critical vesicle size R c for appearance of a new min- 
imum is determined by the condition tp = n or 



ZVi E(R C1 2^/ Rl - Hp = 
8 E(R c ,0) 



(15) 



Substituting the solution of Eq. I|15l) into Eq. © the criti- 
cal surface tension for appearance of the second minimum 
is 



r(l) 
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2/3 



(16) 



^From Eq. |0J we see that a minimum solute concentra- 
tion 



2 (J, 



(i) 



lO^sTA^iWl + ^VK 



(17) 



is necessary to develop the second minimum in AE at r m . 
However, concentration c£2 does not guarantee opening 
of a pore. Even if AE(R,r m ) < 0, the energy barrier 
to pore nucleation can be many fcgT high. Therefore, 
pores with radius less than will quickly re-seal, without 
having a chance to grow. The probability of occurrence of 
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FIG. 1: Energy AE(R,r) necessary to open a pore of radius 
r in a liposome with Ro = 100 nm and ratio 7/(k Ro) — 5 x 
10 _J . The curves are, from top down, R/Ro — 1.0, R/Ro = 
1.0008, R/Ro = 1.001, and R/R = 1.0012 

a sufficiently large thermal fluctuation necessary to open 
a pore with r > is 



P(r) ~ e -/»AB(s,r,) 



(18) 



The waiting time for opening a pore of radius r > r& is, 
therefore, very long unless 



AE{R p ,r b )~k B T 



(19) 



This equation, then, determines the size of a swollen vesi- 
cle 



R p = R Jl + ai 2) /n 



(20) 



which is able to nucleate a growing pore. The membrane 
tension of such a liposome is approximately, 



r(2) 



7T7 

T7r 



(21) 



and the critical pore size is r& ~ 7/c The minimum 
concentration of solute necessary to reach this tension is 



,.(2) 



2 a, 



(2) 



10 3 k B TN A Rp 



(22) 



We note that for membrane parameters given in Table 
1, the membrane tension is — 10 -3 J/m 2 , which is 
very close to the one found to be necessary to rupture a 
mechanically stretched membrane [ll|, [12] • 



III. PORE DYNAMICS 

The rupture dynamics of an osmotically stressed vesi- 
cle proceeds as follows. At t = the vesicle starts 
swelling, its size and internal concentration controlled by 
Eqs. Q an d 0. As it swells, the membrane surface 
tension increase until a — and the energy function 
develops a new minimum. If the barrier height is less 
then k B T, a pore of size is nucleated. From this mo- 
ment the dynamics of the vesicle evolution is controlled 
by the set of Eqs. (QJ, © and Q) . On the other hand, 
if AE(R c ,rb) > k B T, the swelling continues without a 
pore nucleation until Eq. (|19fl is satisfied and a pore of 
radius opens. After a pore is nucleated, the internal 
content of the vesicle begins to leak-out, decreasing the 
membrane tension and leading to an eventual re-sealment 
of the pore. The cycle will be repeated until the in- 
ternal concentration of solute drops bellow c^ in and a 
steady state with j w — is established. In Fig. [21 we 
show the pore radius as a function of time for vesicles 
of three different sizes and initial concentration of solute 
Co = 0.5M. We see that small vesicles are characterized 
by rapid opening and closing of pores, resulting in a peri- 
odic flickering with a characteristic time r/ « 10 _2 s. On 
the other hand, larger vesicles are capable of nucleating 
a long-lived pore. After the long-lived pore has closed, it 
is followed by a sequence of short lived pores, with the 
characteristic life span t/. The life-span of a long-lived 
pore t is show in Fig. [3] For large vesicles, life-span 
scales with the vesicle size as 



(23) 



with v w 2.3 - 2.4, see FigOl 



IV. PHASE DIAGRAM 

To better understand the details of the vesicle evolu- 
tion, it is convenient to separate the membrane dynamics 
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FIG. 2: Radius of a pore as a function of time for vesicles of 
Ro — 200 nm (topmost), 220 nm and 240 nm (bottommost), 
for initial concentration Co = 0.5M. Note that for vesicle 
with Ro = 200 nm the pores are short-lived, while for larger 
vesicles, long-lived pore opens first. 
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FIG. 4: Concentration decay as a function of time. The circles 
are the result of numerical integration of Eqs. Q, © and Q. 
The solid lines are from the analytical expression il lilt . Open 
circles are for vesicle of Ro = lOOnm and filled circles are for 
vesicle of Ro — 300nm. The initial solute concentration is 
c = 0.5M. 
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FIG. 3: Life span of the first-open pore as a function of vesicle 
size Ro for c = 0.5M (rightmost), 1.0M and 5.0M(leftmost). 
Note the appearance of critical vesicle size Rq(co) (sharp 
change in slope of r vs. Ro) which sustains a long-lived pore. 
The peculiar spikes in the life-time of vesicles containing low 
concentration of solute, are an artifact of the way the pores 
are nucleated. On the other hand existence of i?§( c o) is inde- 
pendent of pore nucleation protocol. 



from the concentration dynamics. Since the internal so- 
lute concentration changes very slowly compared to the 
Tf, see Fig. 0] as a first approximation we can take it to 
be constant. 

In this case the vesicle evolution is controlled by 
Eqs. (|TJ and which can be written as 



dR 
~dt 



(24) 



dr 
~db 



G(R,r) 



where 



and 



F(R,r) = 



AirpR 



G(R,r) 



(jw - irr 2 pv) 



1 dE 



rj m £ dr 



(25) 



(26) 



(27) 



The vesicle dynamics is governed by the fixed point 
(r*,R*), determined from dr/dt = = G(R*,r*) = 
and dR/dt = = F(R*,r*). The stability of the fixed 
point is controlled by the eigenvalues Ai and A2 of the 
Jacobian matrix 



J 



dF dF 

OR dr 

dG §G 

dR dr 



(28) 



It is important to keep in mind that the coefficients of the 

Jacobian matrix are real and, therefore, the eigenvalues 

are either real or complex conjugates, Ai = A2. For all 

the parameters that we have investigated the eigenvalues 

are complex conjugates, and the stability is governed by 

Re{\ x ) = Re{\ 2 ) = Re(X). If Re(X) < 0, the fixed point 

is stable and a stationary state with a pore of size r* 

and vesicle of size R* will be established. On the other 

hand if i?e(A) > the fixed point is unstable, and the 

pore will eventually close, see Figs. and El A new 

pore will open when the membrane tension again reaches 
(2) 

the value a c . This process will repeat indefinitely with 
the characteristic time 77. The phase boundary, in the 
concentration-size plane (Ro, c), between the two dynam- 
ical regimes is determined by the condition Re(X) = or 
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FIG. 5: Dynamics of the pore for fixed concentration and 
Ro = 500nm. The top panel is for Co = 0.19M, and the 
bottom panel for Co = 0.18A4". 





t(s) 



FIG. 6: Radius of a pore as a function of time for vesicles of 
Ro = 500 nm for a constant internal concentration of solute 
c = 0.18M (top) and c = 0.19M(bottom). 



equivalently 

which reduces to 
1 d 



TrJ \ R = , 



4irpR 2 OR 



(j w -nr pv) 



1 d 2 E 



r\ m t dr 2 



(29) 



. (30) 



Combining Eqs. J3J and (J5J, the leak-out velocity is 

2ar 

v = — . 

Sirr/R 

Differentiating v with respect to the vesicle size 
dv v I6rn 



(31) 



dR R 3rjA 



(32) 



The change of surface tension with the pore size is given 
by 



da k 
— = -2Trr— 
or Aq 



(33) 



For large concentrations, and R 3> r, the osmotic current 
can be approximated by 



PA. 



so that 



dj w 
OR 



R^ 



At fixed point (R*,r*), 



* 2 * 2PT 
irr pv 



3r/i?o cr* 



(34) 



(35) 



(36) 



where we have approximated R ~ Rq and r* = Tb — 
7/cr*. Eq. (|30|l then reduces to a quartic equation for a* 



27 2 fVmil 



An 



vR 2 



167T T] m £-f 3 K 



. (37) 



Combining Eqs. Q34JI and (|36J) . the phase boundary sep- 
arating the regime of short-lived pores (region II) from 
the long-lived (infinite life-time) pores (region I) is 



Cc{R ) 



PI 



67T7? PRl<J* 



(38) 



R*,r* 



see, Fig0 

For small concentrations, the surface tension does not 
build sufficiently high to cause the membrane rupture 
(region III). The phase boundary between regions II and 
III is denoted by a dashed curve in Fig0 The disconti- 
nuity in slope results from the nucleation barrier passing 
the threshold AE = ksT. Thus the right hand side of 

(2) 

the II-III phase boundary is given by e^-^-Ro), while the 
left hand side c^ n (i?o) 

The role of the concentration dynamics is to make 
the system traverse through the different regions of the 
phase diagram, controlling the time of permanence in 
each regime. 
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FIG. 7: The phase diagram for the constant concentration 
dynamics. In the region I, an open pore is stable. In region 
II, size of a pore oscillates, and in region III pored do not open. 
The dashed line indicates the minimum solute concentration 
needed to open a pore. The discontinuity in its derivative is 
due to the fact that for small Ro the barrier to pore nucleation 
is lower than fc_gT. 



V. DYNAMICS OF SOLUTE LEAK-OUT 

The critical size of a liposome Rq(cq), necessary for 
nucleating a long-lived pore, depends on the initial solute 
concentration. The larger is the solute concentration Co, 
the smaller will be the size of a vesicle which supports a 
long-lived pore. The long life-span of these pores is the 
result of a "wash-out" effect in which the osmotic flux 
is almost completely compensated by the leak-out rate 
of solute through the pore. When solute concentration 
inside the vesicle drops below the critical value, c c {Rq), 
the long-lived pore closes. This value is insensitive to the 
initial solute concentration cq, but depends strongly on 
the vesicle size Rq. 

Combining Eqs.Q and Q we obtain 



4tt 3 dc 
3 dt 



Jv 



Approximating j v 



c(t) = c 



PcA and R 
3Pco 



pRo 



t + 



■ c . 

Ro leads to 
-l 

1 ' 



(39) 



(40) 



Eq. I|4(J|) provides an almost perfect fit of the time de- 
pendence of the internal solute concentration, see FigQ] 
The life span of a long lived pore can be approximated by 
the time it takes for the solute to go from the initial con- 
centration Co to the critical concentration c c (i?o)j below 
which the pore is no longer stable, 



pRo 



3Pc c (Ro) 



1 



Cc(Ro) 



C{) 



(41) 



It is possible to derivate two limits for the critical con- 
centration c c (Rq), see Appendix. Writing 



c c {Ro) = 
for vesicles of radius Rq 



PI 



1 



f(Ro) = 

and for Rq >> R\ 

f(Ro) 



6irT]P R f(R 
i r x = . /aa 

\/3 



(42) 



IT KTj 



31/4 y^" 1 ,s /?( i n i 



-i 2 



Rq ~ R\) 



1 ft 1 / 2 

2 7^ (i? ° 



Ri 



(43) 



(44) 



and therefore the life-span of long-lived pores scales as 



#0' 



R 



2+2/3 



Rq — ^1 
Rq >> i?i . 



(45) 



For the parameters used in this paper i?i = 41.8nm, so 
that Eq. 144(1 is consistent with the numerical findings, 
see Fig. 

The flickering time r/ is approximately the time it 
takes for a vesicle to swell to size R p , needed to induce 
a liposomal rupture. During the swelling, internal con- 
centration of solute changes very little, since R p ~ Rq, 
so that c can be kept constant. Furthermore, for large 
initial solute concentrations, the osmotic current is j w ~ 
A-kPR 2 c 1 and Eq. is easily integrated yielding, 



p{R P ~ Ro) 
Pc 



On the other hand 



(46) 



(47) 



Substituting Eq. gZJ into Eq. (g^ we obtain the expres- 
sion for the flickering time, 



(2) 



p Rq <7, 
2PCK 



(48) 



where a^ 1 is given by Eq. (|2~TJ> . 



VI. CONCLUSIONS 

We have presented a theory for nucleation and growth 
of pores in osmotically stressed liposomal vesicles. The 
model predicts that depending on the internal solute con- 
centration and the liposome size, pores can be either 
short-lived — opening and closing with a characteristic 
time Tf — or long-lived, with their life-time scaling with 
the size of the vesicle. 
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Long lived pores have been observed in red blood cell 
ghosts 0,0]. No theory, up to date, was able to ac- 
count for these long-lived pores. Holes were predicted 
to either grow indefinitely, which would result in ghost 
vesiculation, or to close completely |l5l . Our model pro- 
vides a dynamical mechanism for pore stabilization, con- 
sistent with the experimental observations. However, for 
the specific case of red blood cell ghosts the ratio of -f/n 
must be adjusted to obtain the pore size observed in ex- 
periments. This is not surprising since the real biolog- 
ical cells, unlike liposomes, have a complicated internal 
cytoskcleton, which strongly affects the membrane elas- 
ticity. 

In aqueous solutions the phospholipid membranes ac- 
quire a net negative charge. At physiological concentra- 
tions, 154 mM of NaCl, the Debye length, however, is 
quite short, less then 1 nm and the electrostatic inter- 
actions are strongly screened 0|. We, therefore, do not 
expect that electrostatics will significantly modify the ba- 
sic conclusions of our theory, beyond the renormalization 
of membrane line |l5j and surface tension. However, fur- 
ther, investigations in this direction are necessary and 
will be the subject of future work. 

Finally, up to now we have not taken into account a 
diffusive efflux of solute through an open pore. The char- 
acteristic time for effusion can be estimated as [? ] 



rD 



(49) 



where D is the diffusion constant. Using D w 10~ 9 m 2 /s, 
appropriate for small organic molecules such as sucrose, 
and r — r* sw j/a c w 1 nm, we see that for liposomes 
with R = 200 nm, the time for effusion is r e w 10~ 2 
s. This is comparable to the flicker time r/. Therefore, 
for small vesicles effusion is an important mechanism for 
loss of solute. On the other hand, for large liposomes 
with Rq = 500 nm and above, effusion is only marginally 
relevant. 

This work was supported in part by the Brazilian agen- 
cies CNPq and FAPERGS. 



VII. APPENDIX 



1 (k 

r 2 = - 5 



(53) 



where 



I?! = 



The solutions of a quartic equation like l|51|l can be 
written as 



1 

-1,2 = 5 



-3,4 



V5/± 



2Tx 

Vv 



where y is the real root of the resolvent 

y 3 - 4r 2 y - r 2 = . 

Since the surface tension is non-negative, the physically 
relevant solution for x depends on the sign of Ti and we 
can write 



1 

X= 2 



'2jTi| 

Vv 



\y\ -sgn(ri)^ 



(54) 



For an equation of the form y 3 + ay + b = there is a 
single real root if 



and this real root is 



V = 



1/3 



3) >0 



1/3 



(55) 



This is precisely our case since T 2 < and, therefore, 
Q > 0. It is convenient to rewrite Eg 1551 as 



2 \ 1/3 



Here we present a derivation of the limiting form of the 
phase boundary, Eqs. (|42ll - 144|l . separating the region I 
and II of the phase diagram, Fig. [7| Writing 



7 

a = — x 
no 



Eq. (p?7jl reduces to 



x 4 + T 1 x + T 2 = 



(50) 



(51) 



where 



and 



r 2 \ 3 /4^ 4 



ij){z) = (Vl + z + l) 1/3 - (Vl + z-l)' 
The psi function has the following asymptotic behaviors 



,V3 



with 



* = 4? ( 1 + (52) 



| z -l/3 _ 8 z -4/3 z>>1 



2V3_( Z /2)V3 z<< i 



(56) 
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Putting everything together we have 
r? 

y ■ 



4V 2 



Z » 1 



if 3 z « 1 



In the first case we obtain 



U 



(57) 



and in the second case 



irrl 1 / 3 . 



(59) 



Using the expressions for T% and T2, we obtain the equa- 
(58) tions 
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